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ABSTRACT 

We consider the effects of collective plasma processes on synchrotron emission 
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from highly relativistic electrons. We find, in agreement with Sazonov (1970), 



that strong effects are possible also in the absence of a non-relativistic plasma 
| component, due to the relativistic electrons (and protons) themselves. In con- 

trast with Sazonov, who infers strong effects only in cases where the ratio of 
plasma frequency to cyclotron frequency is much larger than the square of the 
characteristic electron Lorentz factor, v p /vb 3> 7J, we find strong effects also for 
> ; 1 <C v p Ivb "C 7e- The modification of the spectrum is prominent at frequencies 

v < ur* = v v min ^y e , \Jv v )vb j, where vr* generalizes the Razin-Tsytovich 
£sq ■ frequency, u R = 7 e z/ p , to the regime v v jv B <C Je- 

Applying our results to 7-ray burst (GRB) plasmas, we predict a strong mod- 
ification of the radio spectrum on minute time scale following the GRB, at the 
onset of fireball interaction with its surrounding medium, in cases where the ra- 
^p" 1 . tio of the energy carried by the relativistic electrons to the energy carried by the 

magnetic field exceeds ~ 10 5 . Plausible electron distribution functions may lead 

■4-^ 1 

to negative synchrotron reabsorption, i.e to coherent radio emission, which is 
characterized by a low degree of circular polarization. Detection of these effects 



would constrain the fraction of energy in the magnetic field, which is currently 
poorly determined by observations, and, moreover, would provide a novel handle 
on the properties of the environment into which the fireball expands. 

Subject headings: gamma rays: bursts — masers — plasmas — radiation mech- 
anisms: non-thermal — radio continuum: general 



1. Introduction 

According to the model now prevailing (see Piran 2000; Meszaros 2002 for recent re- 
views), Gamma- Ray Bursts (GRB's) originate from the dissipation of the kinetic energy of a 
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relativistically expanding fireball, caused by a cataclysmic collapse of a massive star or by a 
neutron star-neutron star or neutron star-black hole merger event, leading to the acceleration 
of a plasma of electrons and protons to highly relativistic speed. Part of the kinetic energy 
of the this expanding "fireball" is dissipated in "internal" collisions between different parts 
of the inhomogeneous ejecta, resulting in shocks which accelerate particles via the Fermi 
process to ultra-relativistic energies. The non thermal radiation from accelerated electrons 
reproduces well the observed MeV 7-ray spectra. At a later stage of the expansion, the fire- 
ball decelerates due to interaction with its surrounding medium. The relativistic shock wave 
driven into the ambient medium continuously accelerates new electrons of the surrounding 
gas, producing a long-term synchrotron "afterglow" emission. 

Although observations are in general agreement with model predictions (see Kulka- 
rni et al. 2000 for a recent review), the model is incomplete, and there are several open 
issues, which are not properly understood. One is that of the burst progenitor. Sev- 
eral alternative models for the "inner engine" were suggested, such as NS-NS or NS-BH 
merger (Paczyhski 1986; Goodman 1986), and the gravitational collapse of massive stars 
(Woosley 1993; Paczyriski 1998; MacFadyen & Woosley 1999). Unfortunately, at present 
neither GRB nor afterglow observations provide decisive evidence in support of a particular 
model. Nevertheless, the environment may be a clue to the progenitor. Thus, expansion into 
a relatively uniform inter-stellar medium (ISM) with a number density n ~ 1 cm" 3 would be 
a natural consequence of a "merger" scenario, whereas if the progenitor is a collapsing star, 
it is natural to expect a much higher ambient density due to a wind ejected by the star at 
earlier stages of its evolution. The "onset" of fireball interaction with surrounding medium, 
i.e. at the radius where fireball deceleration becomes significant, typically takes place on 
minute time scale (in the observer frame) following the burst, at which stage the density 
of wind plasma is n ~ 10 4 cm" 3 . At present, afterglow observations typically begin several 
hours following the burst and do not allow to directly probe the onset of deceleration. Since 
on a day time scale the fireball expands to a point where the wind density drops to values 
close to that typical for the ISM, present observations do not provide clear discriminants 
between the wind and ISM scenarios (Livio & Waxman 2000). 

A second issue where basic understanding is still lacking is the physics of acceleration of 
electrons to high energies, and the build up of strong magnetic fields by the GRB collisionless 
shock waves (1). The presence of high-energy electrons and strong magnetic fields is implied 
by observations, yet there is presently no theory based on first principles that satisfactorily 
explains electron coupling and magnetic field generation. The ignorance is usually parame- 
terized by two dimensionless parameters, £ e and which stand for the fractions of shock 
internal energy density which are carried by the relativistic electrons and by the magnetic 
field, respectively. Efforts to constrain these parameters using afterglow observations are 
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numerous, and £ e is typically estimated to be close to its equipartition value, i.e. £ e ~ 1/3. 
However, is not well constrained by observations, and its estimated values range from 
£ B ~ 10 _1 (e.g. Waxman 1997a; Wijers & Galama 1998) to £ B ~ 1(T 6 (e.g. Wijers & 
Galama 1998; Chevalier & Li 1999; Galama et al. 1999; Waxman & Loeb 1999). 

Observations strongly suggest that the radiation emitted during the afterglow is syn- 
chrotron radiation. We show here that the modification of the refractive index, by the 
relativistic electrons and protons, may strongly affect the emission at radio wavelengths dur- 
ing the onset of fireball deceleration, if £g <C 1. Moreover, we find that, under plausible 
assumptions on the electron distribution function, the synchrotron reabsorption coefficient 
may become negative, thus leading to coherent emission from the fireball on minute time 
scale following the GRB. These collective plasma processes strongly affect the synchrotron 
spectrum at frequencies < u R * = z/ p y 'v v j v B , where v v fv B ~ a/^/^b is the ratio of plasma 
frequency and electron gyration frequency. 

As mentioned above, currently there is a few hours gap between GRB and afterglow 
observations. However, the (operating) HETE-II satellite and the SWIFT satellite (planned 
to be launched at 2003) may allow observations of early stages of fireball expansion, shortly 
after the GRB, thus providing data on the onset of deceleration. We show below that 
observations of collective plasma effects in the emission from GRB afterglows will serve to 
constrain the value of as well as the parameters of the environment into which the fireball 
expands, thus providing a handle on both progenitor type and shock physics. 

Collective plasma effects on synchrotron emission and reabsorption have been consid- 
ered by many authors (e.g. Ginzburg & Syrovatskii 1965; Zheleznyakov 1967; Yokun 1968; 
Crusius & Schlickeiser 1988). Analyses of effects due to a non- vacuum index of refraction are 
typically limited to the case where the index of refraction is determined by the presence of 
a non-relativistic plasma. These analyses do not apply to the GRB plasma in which we are 
interested, where no "cold" non-relativistic plasma component is present. Sazonov (1969, 
1970, 1973) has studied the case of interest for us, where "cold" plasma is absent, and the 
effects are entirely due to the presence of the relativistic plasma. As we shall show, his 
results are too restrictive, as he infers negative reabsorption only in cases where v p jv B 7^, 
or in cases when the electron distribution function is very anisotropic. We find that these 
constraints might be eased, and that negative reabsorption is possible for isotropic distri- 
bution functions provided that u p ^> ub, thus allowing coherent emission for GRB plasma 
parameters. 

In §2 we derive the dispersion relation and the polarizations and refractive indices of 
normal modes in a highly relativistic, weakly magnetized {y v jvB ^> 1) plasma with isotropic 
electron and proton distribution functions. The frequency range discussed is v v < v <C 
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7e^ p = vr-, where the relativistic plasma may have a strong effect on synchrotron emission. 
In §3 we show that the plasma strongly affects synchrotron radiation at frequencies v < 
v b* — v v\/ v vl u b, which generalizes the Razin-Tsytovich frequency v R to a regime where 
a/ Vp/ v b 7 e - the regime relevant for us. We note, that although v R is commonly quoted 
as the frequency below which plasma effects are strong, it is typically found in numerical 
analyses of non-relativistic plasmas that, for v p /vb ^> 1, synchrotron emission is strongly 
modified only at v < u R * (e.g. Crusius & Schlickeiser 1988). We show in §3 that a more 
careful statement of the qualitative criterion, that leads to the conclusion that plasma effects 
are strong below v R (e.g Rybicki & Lightman 1979 section 8.3), leads directly to the result 
that (for both relativistic and non-relativistic plasma) such effects are important only below 
v R *. 

Our results are applied to GRB plasmas in §4. In §4.1 we briefly describe GRB and 
afterglow phenomenology, give a short review of the fireball model, and derive the plasma pa- 
rameters during the onset of fireball deceleration, considering both expansion into a uniform- 
density ISM (§4.1.1), and expansion into a wind (§4.1.2). Plasma effects are discussed in 
§4.2. The implications of our results are discussed in §5. 

We note, that coherent radio emission from GRB's has recently been discussed by Usov 
& Katz (Usov & Katz 2000). The scenario considered by these authors for GRB production, 
the scattering of ambient medium electrons and protons by a magnetically dominated wind 
(Smolsky & Usov 2000), is different, however, than the scenario we are considering, where ob- 
served radiation is produced (both in GRB and afterglow phases) by the dissipation through 
collisionless shocks of fireball kinetic energy, leading to magnetic field amplification and to 
particle acceleration. The processes we consider are, thus, different than those considered 
by Usov & Katz (2000), and the predicted radio emission is very different. For example, 
while Usov & Katz find power-law spectrum with strong emission at < 1 MHz for strong 
magnetic field, we find strong emission only for weak magnetic fields and over a narrow range 
of frequencies around ~ 0.1 GHz (which is more readily observable). 

Finally , we note that the analysis presented here is restricted to isotropic electron and 
proton distributions, and we show that strong plasma effects on synchrotron radiation occur 
if u p u B . We note however, that anisotropy may lead to further interesting effects, the 
discussion of which is beyond the scope of this work. 
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2. Waves in a weakly-magnetized relativistic plasma with isotropic electron 

and proton distribution functions 

As explained in the next chapter, it is the deviation of the speed at which light propa- 
gates in plasma from c which is responsible for the plasma effects we consider. This implies 
that we must first obtain expressions for the refractive indices rii^iyj) of the transversal elec- 
tromagnetic modes in the plasma. In our derivation of the dispersion relation we restrict the 
discussion to plasmas under the following conditions : (i) Electrons are highly relativistic; 
(ii) Rough energy equipartition between protons and electrons; (iii) Isotropic particle distri- 
bution functions; (iv) Weak magnetization, u p /ub 1. Afterglow observations imply that 
(i) and (ii) hold for GRB plasmas, (iii) restricts the discussion to plasma effects originating 
from the deviation of the refracting indices from their vacuum values. Our discussion is 
restricted to the weak magnetization case since, as shown in the next chapter, strong effects 
of the plasma on synchrotron emission are obtained for this case only. Finally, we restrict 
the discussion to the frequency range affected by collective effects, u B <C v p < v <C r ) e v v . 

The assumption v p /vb ^ 1 allows a perturbative derivation (in vbI v p) of the dispersion 
relation. Moreover, in the frequency range vb <C v v < v <C j e Vp, the deviation of the 
refractive index n = kc/cu, from 1 is much larger than 1/tI- This simplifies the dispersion 
relations obtained below, thus allowing analytic estimates. Below we give order of magnitude 
estimates of the plasma frequency and the difference between the refractive indices of the 
transversal electromagnetic waves in a relativistic plasma. Exact calculations are given in 
the appendix, and are used to verify the simple estimates. 

As we show in §A.l, the refractive indices of the transversal electromagnetic waves in a 
field-free plasma are degenerate, and satisfy approximately the relation 
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where v p is the plasma frequency of a relativistic plasma. This result is exact for monoener- 
getic electron and proton distributions, where 7 e o and 7 P o are the associated Lorentz factors 
of the two species, respectively. Detailed calculations show, however, that the above result 
is a good approximation to the plasma frequency also for more general energy distribu- 
tions. For a power-law energy distribution implied by afterglow observations, 71(7) oc 7~ 2 
for 7 > 7 m i n , replacing 70 with 7 m j n in equation (1) gives the plasma frequency with accu- 
racy of 1%. Note, that equation (1) is a natural generalization to the relativistic regime of 
the familiar expression for the plasma frequency, obtained by replacing the particle mass m 
with 7m. Since we assume that the energy of the electrons and protons is approximately 
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equipartitioned, the contributions of the two species to v v are comparable. 

Once an external magnetic field B is introduced, the plasma cannot be regarded as 
isotropic any longer, and the degeneracy of the two refractive indices is removed. For fre- 
quencies v <C ub the effect of the magnetized plasma on the propagating radiation is small, 
and the electric field can be assumed to be approximately transversal to the direction of 
propagation, i.e. E _L k. In the appendix (§A.2) we derive the dispersion relation for the 
waves of transversal electric field. In order to obtain estimates of the refractive indices of the 
transversal EM modes, we solve the Vlasov equation of a relativistic plasma with isotropic 
electron and proton distribution functions. The assumption of weak magnetization allows a 
perturbative expansion of the equation, where the perturbations in the particle distribution 
functions are assumed to be linear in the external magnetic field, (i.e., the expansion param- 
eter is u B /u p ). This leads to a significant simplification of the formalism. The dispersion 
relation for the components of the transversal electric fields is given by : 



where u = 2tiv is the radian frequency of the wave, k is the wave-number, and is the 
2x2 dielectric tensor. The refractive indices n li2 (^) = ck lf2 /uj of the plasma are the square 
roots of the eigenvalues of this equation. Detailed numerical calculations showed that the 
values obtained for the components of the dielectric tensor and the resultant refractive in- 
dices are not sensitive to the exact shape of the electron and proton energy distribution 
functions. Specifically, we showed that simple expressions obtained for monoenergetic distri- 
bution functions (which are the estimates used below) are good approximations of the more 
cumbersome expressions obtained for power law distributions. 

The deviation from degeneracy of the two refractive indices is the consequence of the 
non-diagonal terms of the dielectric tensor, which are proportional to cos0, where is the 
angle between the wave- vector k and the external magnetic field B [eqs. (A5a) and (A9) ]. 
Hence the deviation is largest for radiation propagating along the magnetic field (0 = 0), and 
vanishes when the direction of propagation is perpendicular to the magnetic field (0 = 7r/2). 
Equations (A6) and (A10), describing the diagonal and off-diagonal elements of the dielectric 
tensor, can be used to estimate the difference between the two refractive indices. For the 
frequencies in which we are interested 1 — n is larger, or at least comparable to I/7I p . 
This simplifies the algebra immensely. Since we assume that the widths of the electron and 
proton distribution functions are close, and that energy density in the protons and electrons 
is approximately equipartitioned, the contributions of the two species to An = n\ — 112 are 
comparable. Hence, an order of magnitude estimate of An is : 
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up to a factor of order 1. 

The result An ( f )=n /2 = is a consequence of the fact that our derivation keeps only 
terms which are linear in the magnetic field, i.e. 1 st order terms in u B /u p . If this assumption 
is eased, and the refractive indices of the ordinary and extra-ordinary modes are used to 
estimate An0 =7r / 2 , we obtain An^, =7r / 2 = \{y B jv) 2 [l — v 2 /2v 2 ~\. As will be shown in the 
next chapter, the frequency at which the plasma has a significant effect on the synchrotron 
emission is v B * — v p\Z v pl v b- This implies that for plasma parameters of interest to us 
{y v jv B ^ 10 3 ) the logarithmic factor in equation (3) is of order 10. Hence we replace 
equation (3) with an expression which takes into account the finite An at perpendicular 
propagation : 



An(0) 
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(4) 



otherwise. 



We solved the dispersion relation [eqs. (A6) and (A10)] numerically, for two values of 7 e 
(920 and 4.6 x 10 4 ), and for three values of u B /u p (10 -1 , 10~ 2 and 10~ 3 ). These values of 7 e 
and v B jv v are representative of typical values for GRB afterglow plasmas (see §4). In these 
calculations we assumed a propagation angle = ir/4. The results agree with the estimate 
(4) up to a factor of ~ 2. The numerical calculation also confirmed that for parameters in 
the ranges relevant for GRB afterglows, the discrepancy between the refractive indices of the 
transversal modes is small in comparison to the deviation of either of them from unity, i.e., 
\ni - n 2 \ < (1 - ni )2 ). 

When the radiation propagates parallel to the magnetic field, the squared refractive 
indices of the two normal modes may be described by the well-known expressions 

v 2 v 2 

, _ 1 P e PP 

£r — 1 — 



v{v - v Be ) v{y + v Bp ) 
2 2 

V V 
V\V + "Be) V[y - V Bv ) 

where v Be ,v Bv are the gyration frequencies of electrons and protons in the magnetic field, 
and v pe , Vpp are the relativistic plasma frequencies of the two species, respectively. It is 
instructive to compare our estimate in equation (4) to the difference n R — obtained 
from equation (5). When exact equipartition is assumed, and both the protons and the 
electrons are mono-energetically distributed, u pe = v vv = v v and u Be = u Bp = v B , whence 
e B = €l = 1 — 2v 2 j{y 2 — u B ). This degeneracy is removed if equipartition is not exactly 
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satisfied, or if the widths of the distribution functions of the two species are different. Since 
both effects are expected to be important in the case of our interest, we will use equation 
(4) as an estimate for An. 

Our calculations imply that non-diagonal elements of the dielectric tensor are smaller 
than its diagonal elements by a factor of the order of v B jv <C 1. The latter, however, are 
independent of the magnetic field, and thus are identical to the values we obtained for a 
field-free plasma [see eq. (A5a)]. As will be discussed in §3, the difference between the 
refractive indices of the transversal modes is important for determining the polarization of 
the modes and thus for the calculation of synchrotron self-absorption. However, our estimate 
of the frequency at which plasma effects on synchrotron emission become significant is not 
sensitive to An. Therefore, as far as estimating the frequency u R * = v p ^Jv v jvB is concerned, 
we will use the approximation given in equation (1). 

The detailed expressions we obtained (in the appendix) for the elements of the 2x2 
dielectric tensor e^, show that en = €22 and 612 = —621 [see eqs. (A6) and (A10)]. The 
assumption of transversality thus implies that the normal modes are always circularly po- 
larized (clockwise and anti-clockwise). To check the consistency of our assumption of quasi- 
transversality, and to obtain the (elliptical) polarization of the two quasi-transversal modes, 
we calculated the eigenmodes of the full 3x3 dielectric tensor. Following the method outlined 
in §A.l [and described in the paragraph preceding eq. (2)], it is straightforward (yet lengthy) 
to obtain expressions for the other elements of with i,j = 1,2,3, and without assuming 
transversality. It is possible to show that these extra elements contribute an additional term 
to An, which is smaller than the value obtained from the transversal calculation by a factor 
of ~ 0.1(ln |i/ s /4i/p|)" 4 z/ p /z/ B ~ 1CT 2 at small and mildly oblique propagation angles, and 
becomes comparable to | (v B /is 2 ) [l — (z/ 2 /2z/ 2 )], as \n/2 — <p\ approaches 0.1 v B vjv 2 p . Nu- 
merical calculation of the plasma normal modes confirmed that for plasma parameters of 
interest to us (7 e = 10 3 — 10 4 , u p /u B = 10 2 — 10 3 ) the normal modes are indeed transverse 
to the direction of propagation, and are left- and right- circularly polarized. 



3. Effects on synchrotron emission 

The spectral characteristics of synchrotron emission by a single relativistic electron are 
determined by the beaming of the radiation emitted by the electron into a narrow cone about 
the instantaneous direction of the electron's motion. When the electron is in a vacuum, the 
opening angle of this cone depends on the electron's Lorentz factor and the frequency of the 
radiation. Three regimes should be considered: when the frequency of interest is higher than 
the characteristic synchrotron frequency v c the cone's opening angle is A9 ~ (z/ c /z/) 1 / 2 /7, and 
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when v ~ v c , one obtains the familiar result A6 — I/7; If, however, the frequency of interest 
is much smaller than v c , one obtains A9 ~ (y c / z/) 1 / 3 /j. Next we consider an electron which is 
embedded in a dielectric medium with a refractive index 1. Since the speed of light is now 

1 /2 

c/n, the opening angle of the cone is A9 ~ [(z/ c /z/)(l/7 2 ) + 1 — n 2 (u)] for v > u c , A9 ~ 
[I/7 2 + 1 — n 2 {y)] 1 ^ 2 for frequencies v ~ z/ c , and A9 ~ [(z/ c /V) 2 / 3 (l/7 2 ) + 1 — n 2 (z/)]^ 2 for 
frequencies much smaller than v c . When the dielectric medium is a plasma, we use 1 — n 2 = 
Vp/v 2 ■ A strong deviation from the vacuum behavior occurs when the second term iy v jv) 2 
is comparable to the first term in the square brackets. Using v c = 7 3 z/_b, where ub is the 
gyration frequency of the electron, we obtain an expression for the frequency at which the 
plasma collective effects greatly influence the synchrotron emission : 



v R * ~ u p min |7, J — j . (6) 

The Razin-Tsytovich frequency vr = Vp'-f is usually given as an estimate to the frequency 
below which the plasma strongly affects synchrotron radiation. Equation (6) shows that this 
estimate is not generally valid. We are interested in the regime sj Vp/vs 7 e , hence for 
plasma parameters prevailing in GRB afterglows the frequency at which the plasma greatly 
affects the synchrotron radiation is vr* = v v \/ v v j vb- It is simple to show that the definition 
of the plasma frequency [eq. (1)] implies that = Ce/^B, where £ = \n(^ emax /^ emin ). 

Typically I ~ 4. Thus 

1^ V4 



VR * = Vv H~B = Vv \nT B ) ■ (7) 

An immediate consequence of equation (7) is that for the effect to be apparent we must 
have Up 3> v B (or, equivalently, £ e 3> 2££b)- Sazonov (1970) also considered the effects 
of a relativistic plasma on synchrotron emission. However, he considered only a situation 
in which v ~ v c 7^, which led him to infer that negative synchrotron reabsorption is 
possible only for v v jv B ^> 7 2 , or equivalently, if ^ e /Cs ^ it- The result (7) shows that this 
condition is far too restrictive. However, as will be discussed in §4.2, for the effects to be 
within the detection range of current radio telescopes, we must have £ e /£e ^ 10 5 . 



3.1. The Razin suppression 

All the information regarding the response of the plasma to electromagnetic waves (e.g. 
modes propagating through the plasma, loss and gain, etc.) is contained in the dielectric 
tensor. Specifically, the strong effect of plasma on synchrotron radiation is a manifestation 
of the deviation of the refractive index from 1, that is the non-trivial (i.e. "non- vacuum" ) 
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structure of the dielectric tensor. The expression derived in the appendix for the dielectric 
tensor is not complete, i.e. it does not include self-absorption, since we treat the magnetic 
field perturbatively. Unfortunately, obtaining a complete expression for the dielectric tensor 
requires solving a kinetic equation for the plasma, a process which may be avoided by 
introducing the method of Einstein coefficients. The application of this method requires that 
over the length-scale of the wave phenomenon (i.e. the wavelength), the absorption does not 
change the wave characteristics considerably. Indeed, since we do not obtain the dielectric 
tensor explicitly, we do not know a priori the EM modes which can propagate through 
the plasma. We must assume, however, that whatever these modes are, the modifications 
introduced to the dielectric tensor due to the presence of the plasma are dominated by the 
deviation of the refractive indices from unity, and not by the absorption, that is : 



where A = cjv is the wavelength, and / denotes the mode. When condition (8) is not 
satisfied one must use a kinetic-equation approach, to self-consistently derive a dispersion 
relation which incorporates the self-absorption. 

Since the refractive index depends on the radiation mode under consideration, we now 
have to determine the modes relevant to the frequency range we are interested in. Motivated 
by the insight that in "vacuum" absorption is present, and the two transversal modes share 
the same refractive index (i.e., 1), whereas, on the other hand, the result obtained above 
(see §2), shows that when absorption is neglected one obtains a finite An® [see eq. (4)], 
we conclude that the question of determining the relevant modes is resolved by considering 
whether the following inequality holds : 



where a v is the synchrotron self- absorption coefficient, and A the wavelength of radiation. 
When condition (9) is satisfied, the normal waves are linearly polarized along and perpendic- 
ular to the projection of the magnetic field on the plane of observation (Ginzburg 1989). We 
denote these polarizations by _L and ||, respectively. The power emitted by a single electron 
is then given by 




(8) 




(9) 
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where B is the strength of the magnetic field, and x the pitch angle. 

If, however, condition (9) is not satisfied, one cannot neglect the difference between the 
refractive indices of the normal modes. We are thus required to consider the two normal 
modes of the plasma, discussed previously in §2. The normal modes are (quasi-) transversal, 
and are left- and right- circularly polarized. It can be shown (Ginzburg 1989) that half the 
total power is "converted" into each circularly-polarized normal wave (if we take into account 
terms up to order I/7). In this case, the power emitted by a single electron is 

P^(,.T) = ^gf^ [l + 7 2 (l-< 2 )]- 1/2 ^ f Wfe, (11.) 

ZeBsinx 2 r, . 2/ n 2 x -1-3/2 ( s 

*c = -^ — 7 [l+7(l-n li2 )] • (Hb) 

The polarization of synchrotron emission in plasma will be further discussed in §3.3. 

If the frequencies of interest are higher than vr* , the emission does not differ greatly from 
the emission in vacuum. However, for frequencies v < vr* the emitted power is dramatically 
suppressed for both polarization regimes. This suppression is usually called the Razin effect. 



3.2. Synchrotron self-absorption and the possibility of negative reabsorption 

Our calculation of the synchrotron self-absorption coefficient a u implements the Einstein 
relations between emission and absorption coefficients. These relations relate the emission 
and absorption coefficients of photons having a specific polarization state, or, to put it in 
terms of the previous chapter - a specific mode propagating through the plasma, which we 
shall denote by I. Hence also the expression for the self-absorption coefficient is for a specific 
mode (Ginzburg 1989) : 

where P( l \v, E) is the power of /—polarized photons emitted by an electron with energy E. 

Negative contributions to reabsorption come only from regions where the electron dis- 
tribution function grows faster than E 2 . Indeed, if the distribution function has regions 
which are "steep" enough, the self-absorption coefficient becomes negative at low frequen- 
cies [y < vr*). This corresponds to stimulated emission from the plasma, and radiation is 
coherently amplified as it propagates. For this reason the effect is sometimes known as a 
maser effect. 
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In order to estimate the "amplitude" of the negative self-absorption, we take the electron 
distribution to be mono-energetic. Substituting n e (E) = n e 5(E — E Q ) in equation (12), and 
using the fact that we are interested in frequencies u <C 7 e f p , it can be shown that 

a u ~2x 10- 2 (u p u B /c^f e u)^(Z(u)) , (13) 

where $(Z) = 2Z f™ K 5/3 (y)dy -2ZK 5/3 (Z), and Z{v) = |^/z/ B z/ 2 . The function $(Z(z/)) 
has a global minimum at approximately u R * = v v ^/ v v j vb, and its minimal value is —0.24. 
Substituting v ~ ur*, we obtain an estimate of the minimal value of the synchrotron self- 
absorption : 



mina 




up to a factor of order unity. 



3.3. Polarization of synchrotron radiation in plasma 

We can now use the estimate of the absorption coefficient obtained above to check the 
self-consistency of our calculation. Namely, we shall show that in the frequency range of 
interest, the deviation from unity of the refractive indices dominates the absorption, thus 
facilitating the validity of implementation of the Einstein coefficient method, as required 
by the consistency condition (8). In §A.l we showed that the refractive indices may be 
approximated by n ~ 1 — \{y v jv) 2 . Hence, for 5— function distributions, equation (8) may 
be re-written as 10 2 y/%Up^ 3 / u B ^ 2 ^> v. Since we are interested in frequencies of order of 
v R * ~ Vp^/vp/vB-, this condition is casted to 10 2 ^/t^z/ p /z/b 1, which is always satisfied for 
typical afterglow parameters (7 e ~ 10 3 , v v jv B ~ 10 3 ). 

Next we determine the polarization of the synchrotron modes in the frequency range 
of interest for us. We use the estimate (4) for An, and substitute equation (14) for the rhs 
of the inequality (9) with negative reabsorption. After some algebra we find that condition 
(9) is satisfied only for radiation propagating almost perpendicular to the magnetic field, 
\ir/2 — <f)\ < 0.1(vb/v p ) 1 ^ 2 given that v B jvp 0.02'je 1 ^ 2 - For typical plasma parameters 
(7 e ~ 10 3 — 10 4 , Up/ u B ~ 10 3 ) this inequality is not satisfied, and we conclude that the plasma 
is not isotropic, hence the normal waves at frequencies v ~ u R * are circularly polarized. 



4. Application to GRBs 



This section is dedicated to the study of the application of the results obtained above 
to the plasma conditions prevailing in GRB afterglows. We briefly review the fireball model, 
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and obtain the physical parameters of the plasma during the onset of fireball deceleration, 
considering both expansion into a uniform density ISM (§4.1.1), and into a wind (§4.1.2). 
These parameters enable order of magnitude estimates of relevant frequencies (§4.2.1 and 
§4.2.2). We then make a short excursion, referring to the issue of the origin of negative 
reabsorption in relation to the shape of the electron distribution function at low energies 
(§4.2.3), and conclude the section with a numerical calculation of the resultant spectrum, 
for various electron distribution functions (§4.2.4). 



4.1. Fireball plasma parameters during the onset of deceleration 

According to the fireball paradigm (Paczyhski 1986; Goodman 1986), an energetic ex- 
plosion (E ~ 10 51 — 10 53 ergs) drives a relativistic blast wave into an ambient gas (the 
"forward shock"). The expanding shell of accelerated ambient gas approaches gradually a 
self-similar behavior (Blandford & McKee 1976). During the short transition period before 
self-similarity is established, the interaction between the ejecta and the surrounding medium 
drives a relativistic shock into the ejecta (the "reverse shock") (Meszaros & Rees 1997), and 
heats it. The transition to self-similarity occurs on a time scale comparable to the time it 
takes the reverse shock to cross the ejecta (e.g. Waxman & Draine 2000). The observed radi- 
ation is the consequence of synchrotron emission by relativistic shock-accelerated electrons, 
which gyrate in the magnetic fields generated by the shocks. 

Recent data support GRB models where the outflow is a jet (rather than a sphere), with 
an opening angle 9j et ~ 0.1 (Waxman, Kulkarni & Frail 1998; Fruchter et al. 1999; Stanek et 
al. 1999; Harrison et al. 1999; see Frail et al. 2001 for update analysis). The dynamics and 
resulting light curves of such models differ from those of isotropic expansion models after 
the jet Lorentz factor decreases below \jQj et . Typically this happens several hours after 
the main GRB. We, however, are interested in very early stages of the expansion (typically 
~ 10 s after the main GRB), and so the analysis we present below, although formulated in 
terms of an isotropic model, is valid also for a jetted scenario. 



4-1.1. Expansion into a uniform- density ISM 

Self similarity is established once the reverse shock crosses the ejecta. It has been 
shown (Waxman & Draine 2000) that at this time both the shocked ISM and the heated 
ejecta propagate with a Lorentz factor which is close to that given by the Blandford-Mckee 
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self-similar solution : 

where T = lOTi s is the observed burst duration, which is typically of the order of 10 seconds, 
E = 10 52 E 52 ergs is the explosion energy, and n = Iuq cm -3 is the density of the ambient 
gas. Accordingly, the electron number density behind the forward shock is 

n KF) = 4T (F) n = 73 5£V8 Ti -3/8 n 7/8 cm _3 ^ 

(where the prime denotes that this is a quantity measured in the frame comoving with the 
plasma). 

Let 7p and 7p be the Lorentz factors associated with the thermal motion of the 
protons accelerated by the forward and reverse shocks, respectively. Then it can be shown 
(Waxman & Draine 2000) that 7^ ~ 7p R - ) r 2 /r i , where Tj ~ 300 is the Lorentz factor of 
the ejecta prior to its deceleration by the ambient gas (and after the production of the main 
GRB). This result holds for both relativistic and non-relativistic reverse shocks. Since the 
plasmas behind the forward and reverse shocks are separated by a contact discontinuity, the 
energy densities in the two plasmas are similar. Consequently we obtain 

nf ) ~ (T 2 /T t )nf > = 8.27 x 10 4 £ 5 3 2 /8 Tr 9/8 nf cm" 3 . (17) 

Lacking a fundamental theory of Fermi acceleration and formation of magnetic fields by 
shocks, it is customary to parameterize the fractions of the energy carried by the magnetic 
field and the electrons by two dimensionless parameters, and £ e , respectively. We assume 
that the values of these parameters are similar in the plasmas accelerated by the forward 
and reverse shocks, as we expect them to be associated with "micro - physics" processes in 
the plasma. Since energy densities behind both shocks are similar, we then have 

B (R) „ B (F) = 0.07 (^_^ V2 SV^/S^/S q (lg) 

where we use the normalization (£s/10~ 6 ) since we expect strong collective plasma effects 
for small values of £g (see §1). 

Relativistic shock waves are assumed to accelerate protons and electrons to high energies, 
giving power-law distribution functions : 



7^(7) = if 7 P , for 7 emin < 7 < ^ e max 



(19) 
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Spectral indices p > 2 were observed in cosmic-rays spectrum and supernovae radio emission, 
observations which were later explained by theoretical work (Blandford & Eichler 1987). 
Numeric and analytic calculations of particle acceleration via the first-order Fermi mech- 
anism in relativistic shocks yield spectral indices p ~ 2.2 for highly relativistic shocks 
(Bednarz & Ostrowski 1998). This result is in agreement with the value of p inferred from 
GRB afterglow observations (Waxman 1997a, 1997b). The maximum Lorentz factor ^emax 
is determined by requiring that the most energetic electrons lose energy through synchrotron 
emission slower than they gain energy due to acceleration by the shock. Normalization then 
implies 

7e2 = = 4.22 x 10V (^) E^T-^nT , (20) 

where t = e max H emin)- Typically t ~ 4. Following the discussion leading to equation 
(16), we have 7,^7^/^. 



4-1.2. Expansion into a wind 

We examine here the simplest model of a "wind" , in which the star loses mass at a con- 
stant rate M during an epoch prior to the explosion. Material is ejected radially at a constant 
speed v, hence producing a non- homogeneous ambient gas with n oc r~ 2 . Using typical es- 
timates for the mass- loss rate M = lO~ 5 M /year and for the wind velocity v = 10 3 km s~ x 
(Chevalier & Li 1999), one obtains n = M / 4iim p vr~ 2 ~ 3 x 10 35 ^M_ 5 /t> 3 j r~ 2 cm -3 ; here 

we used the notations M_ 5 = M/(lO 5 M /year), and t> 3 = t>/(10 3 km s _1 ). Once the dy- 
namics reaches the self-similar regime, the Lorentz factor and the radius are related through 
E = (16-7T /9)m p c 2 r 3 n(r)r 2 B _ M = (4Mc 2 /9v)T 2 B _ M r. The physical parameters of the plasmas 
heated by the forward and reverse shocks are obtained following the reasoning employed in 
§4.1.1, and so, without further ado: 

n'P ~ 1.09 x l0 7 E- 2 3/4 (M_ 5 /v 3 ) 7 / 4 T- 5/4 cm" 3 , (21) 
nf) ~ 6.56 x 10 7 E-V\M_ 5 /v 3 ) 5 '*T- 7/A rrl 5 cm" 3 , (22) 

B m ^ ^)^4.1 8 (^y /2 ^ /4 (M^ 3 ) 3/4 Tr 3/4 G, (23) 
7e <£ - 9.78 x 10V (^) ^ 2 /4 (M- 5 M)- 1/4 ^ 1/4 , (24) 
7.2. - 1-61 x 10V (J^j ^(M^) 1 ^ . (25) 
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4.2. Razin cutoff and synchrotron maser 

We shall now use the lessons of the previous three chapters to estimate u p , ub, v c (the 
synchrotron characteristic frequency) and vr* for the four cases under study, i.e. the ejecta 
heated by the reverse shock, and the ambient medium heated by the forward shock, in the 
two fireball expansion scenarios. 



4-2.1. Expansion into a uniform- density ISM 

The typical Lorentz factor of the ISM accelerated by the forward shock is very high, 
lemin — 4-2 x 10 4 , and the number density is n'i F ' ~ 735 cm -3 . Under the assumption of a 
small £b, the magnetic field is B ~ 0.07 G. The values just stated correspond to the epoch 
of transition to self-similar behavior. Using the results of the last section, and boosting to 
the observer frame (multiplying by T), we obtain : 

'1/2 



i/f> ~ 2.2xl0 5 4 /2 ^" EH 8 T- 3/8 n 3 /8 Hz 
VB " 8 ^(^)" 1 (^)^^^ /Sno " 1/BHz 



v<P - 9.8 x lO^if (||) 2 (J^j V2 EU%-^ Hz 

i/g? - 3.5xl0 6 4 /4 (||)" /4 (^)" 1/4 ^Tr 3 ^ (26) 

Since 7^ in , -B and fie F ' are proportional to T oc r -3 / 2 , we find that (in the observer's frame) 
v p , vb and vr* oc r~ 3//2 , whereas z/ c oc r~ 6 . 

The ejecta heated by the reverse shock is characterized by a lower Lorentz factor 7 e „ li „ — 
375, but the number density is much higher than behind the forward shock, nr' ^ 8.3 x 
10 4 cm -3 . In the observer's frame : 

*f> ^ 2.5 x we {^) 112 Etl%-v\l' 8 r-l 5 m 

/ t \ - 1 / t \ V2 



^ ^ 7.7X10 1 V(||) ( 



2 / t \ V 2 

1 7ln /2 r 2 9 = Hz 



10-6 / 11 ' ''2.5 

-1/4 / c n -1/4 



i/g> - 3.9xl0^ 4 (i|)" (^)~ ^rr 9 ^^ Hz . (27) 
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4-2.2. Expansion into a wind 

While the reverse shock crosses the ejecta, the heated wind particles are accelerated by 
the forward shock to a typical Lorentz factor 7 e ^^ — 9.8 x 10 3 . At that time the number 
density is n'i F ' ~ 1.1 x 10 7 cm -3 , and the magnetic field is B ~ 4.2 G. Then : 



vjP * 1.3 x m\' 2 (j^y (Af_ 5 M) 3 W /4 Hz 

v<P - 5.1 x 10% (||) " (J^J V2 {M^/v^E-^T-^ Hz 

i/f> - 7.1xl0^ 2 ^ 2 ^y /2 ^rr 3/2 Hz 



(F) ~ 2.0 x 10% /4 (Jl) (M- 5 M) 3/4 ^ 2 1/4 Tr 3/4 Hz . (28) 



V R* 



Since during the self-similar stage n' e scales as T 5 , B tx T 3 and 7 em m oc T, and T oc r -1 / 2 , 
we obtain the following scaling relations for later times : u p , ub and u R * oc r -3 / 2 , z/ c oc r~ 3 . 
However, for earlier times the forward shock expands at a uniform Lorentz factor, and so 
lemin oc T = const., n' e oc Yn oc r~ 2 , and -B oc v^w 7 oc Vr 2 n oc r~ l . Hence we find that 
VpyVs, v c and v R * all scale as r _1 . 

And finally, the electrons in the ejecta heated by the reverse shock when the fireball 
expands into a wind have a Lorentz factor of - 1.6 x 10 3 , and number density n'J^ ~ 
6.6 x 10 7 cm -3 . Hence : 



4 R) 



■ 7.7 x 10 7 £f (J^j ' (M_ 5 /v 3 )^%~ 5/4 T- 2 \ Hz 

<?> - 3.1 x 10% (|£f (^y^ (M^M) 1 / 4 ^ 5 / 4 !- Hz 

W - 2.0 x 10 15 C 2 (||) 2 (^) V2 (M_ 5 M)^ 2 1/2 Tr 1/2 r 2 , 5 Hz 

g> - 1.2 x 10YJ 4 (i|)" 1/4 (^) _V \M_ 5 M) 1/4 ^Tr 5/4 r- (29) 



4-2.3. The electron distribution function at low energies 



As stated in §3.2, a necessary condition for negative reabsorption is a region which grows 
faster than E 2 in the electron distribution function. We give here a plausible mechanism 
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which may be responsible for the existence of such a region. Relativistic shock waves are 
thought to accelerate electrons in such a way that the electron distribution function has a 
power-law tail extending to high energies. However, there is no satisfactory theory which 
predicts the shape of the distribution function at low energies. We thus assume that at 
low energies the distribution function has the simplest possible behavior, and electrons dis- 
tribute according to the volume they occupy in phase-space, i.e. n e (E) oc E 2 . By emitting 
synchrotron radiation, the electrons gradually lose their energy, and accumulate at lower 
energies, thus leading to a distribution function which grows faster than E 2 at low energies. 
The excess of electrons above the E 2 power law is sensitively dependent on the details of the 
distribution function injected by the shock. Consider, for example, an injected distribution 
function composed of two pure power laws [i.e., n e (E) oc E 2 if E < Eq, n e (E) oc E~ p if 
E > E ]. This distribution function has a discontinuous derivative at its peak. Synchrotron 
cooling thus results in an excess of "cooled" electrons just below E , leading to a strong 
maser effect. This issue will be treated further in §4.2.4. 

It must be borne in mind, though, that electrons which are "injected" to the hot plasma 
by the shock at later times have less time to cool. This leads, in principal, to a significant 
complication in the distribution function of the shock-accelerated electrons, since the further 
away from the shock front we look, the "colder" the distribution is, and in general, we cannot 
consider the distribution function of the shocked plasma to be homogeneous. In particular, in 
regions closest to the shock front, where the electrons had very little time to cool, compared 
to the dynamical time, the reabsorption is positive, and thus may obscure coherent emission 
from regions further away from the shock front. Nevertheless, we have shown numerically 
that these regions contribute an optical depth of order of few, at most, for typical afterglow 
parameters; hence the effect of inhomogeneity does not change qualitatively the phenomenon 
of coherent emission by the relativistic plasma. Furthermore, if one assumes some turbulent 
mechanism which acts on the dynamical time scale of the system, and "mixes" electron 
populations which were injected by the shock at different times, one may disregard that 
complication, and treat the "averaged distribution" function. 

It is interesting to notice that a major effect exists even when the low-energy electrons 
are distributed as n e (E) oc E 2 ! Indeed, there is no negative reabsorption in this case, since 
a v < requires a region steeper than E 2 in the distribution function. Nevertheless, the 
emissivity at v ~ ur* is dominated by the low-energy electrons (since j u of electrons with 
7e > lemin is Razin-suppressed at higher frequencies). On the other hand, the absorption 
coefficient, albeit always positive, is dominated by the high energy electrons, and so starts 
decreasing at frequencies higher than u R * . Altogether, we obtain a high peak in the intensity 
for t u 1, and a strong suppression when r v <C 1, where r v „ is the optical depth at vr*. 
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4-2.4- Results of detailed calculations 

We calculated the emitted intensity for a plasma with 'jemin = 10 3 — 10 4 , v v jvB = 
10 2 — 10 3 . Following the discussion in §3.3, the normal waves propagating in the plasma are 
circularly polarized. The intensity of the emitted radiation is given by 

1 — e -7 "^ 1 - 2 ) 

4(i,2) = Mi,2)&' , (30) 

7V,(i,2) 

where j u is the specific emissivity, and t v the optical depth. The width (as measured in the 
comoving frame) of the emitting medium along the line of sight A' is typically ~ 10 13 cm in 
both expansion scenarios. 

Using a shock injected electron distribution function consisting of two pure power laws 
with typical afterglow parameters, we verified that the cooling scheme suggested in §4.2.3 
results in a maser effect at frequencies v < v R *. More realistic distribution functions are 
expected to have a smooth transition region between the two power laws. As an example, we 
considered the following shock injected electron distribution function (following Gruzinov & 
Waxman 1999): 

n' e (z) txz^l + az-^y 1 , (31) 

where z = J fremiti- The constant a and the overall proportionality factor are chosen such 
that the distribution is adequately normalized. This distribution function, behaves asymp- 
totically as 7 2 at low energies and as 7~ p at high energies, with p = 2.4. We applied 
the cooling scheme described in §4.2.3 to two sets of plasma parameters : the one, with 
v p jv B = 250,7 emin = 4 x 10 4 ,rig = 10 3 cm -3 and 'j coo i/le m in - 2 x 10 3 , and the other, with 
v p jv B = 250,7 ernira = 10 4 ,n' e = 10 7 cm" 3 and j cool ~ I0j em in- Here 7 coo/ = 6nm e c/o- T B 2 t' dyn 
is the maximal electron Lorentz factor allowing significant synchrotron losses within the dy- 
namic time scale of the system, t' d . The first parameter set corresponds approximately to 
the conditions prevailing in the ISM heated by the forward shock, when the fireball expands 
into a uniform-density ISM, while the other corresponds approximately to the physical con- 
ditions in the plasma accelerated by the forward shock, when the fireball expands into a 
wind. The results of the numerical calculations are displayed in Figure 1. When the first 
parameter set is used (Figure 1 Left), there is no negative reabsorption, since the cooling is 
not efficient enough. An increase of ~ 1.5 orders of magnitude in the emitted intensity, com- 
pared to the emission in "vacuum", is apparent, as explained in §4.2.3. On the other hand, 
the second parameter set, which corresponds to a much lower value of r y C ooi/lemin, (Figure 
1 Right) clearly shows the maser effect at frequencies v < z/r*, reflecting the negative self 
absorption coefficient due to the efficient synchrotron cooling. 
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5. 



Discussion 



We have derived a dispersion relation for transversal electromagnetic waves in a weakly 
magnetized relativistic plasma, assuming the electron and proton distribution functions are 
isotropic. The frequency range at which we were interested was u p < v <C 7 e z/ p . Treating the 
external weak magnetic field as a perturbation, we have shown that at this frequency range 
the normal modes are circularly polarized, having refractive indices which can be approxi- 

1 /9 

mated by the familiar expression n\ 2 ~ 1 — (z//V p ) 2 , where v v = \{Vp e R ) 2 he + iyp P R ) 2 / l-p\ 
We have shown that this result is valid whether the electron distribution is assumed to be 
a 5— function or a power-law. We have also calculated the difference between the two re- 
fractive indices [see eq. (4)], obtaining |Ara| ~ (z/|z/b/V 3 ) cos0 In (z/ 2 /4z/ 2 ), for radiation 
propagating at angle <fi < \n/2 — 0.1z/bz//V 2 | with respect to the external magnetic field, and 
|An| ~ v^jlv 2 for radiation propagating at larger angles (i.e., almost perpendicular to the 
field direction). 

We next came to consider the effects of the relativistic plasma on synchrotron emission. 
We have derived an estimate of the frequency at which the plasma effect on the emission 
becomes significant : 



I V u b J 

This result generalizes the familiar Razin frequency vr = 7 e z/ p to a regime where v p /vb <C 
7g, which is the relevant regime for typical plasma parameters at GRB afterglows; hence 



we expect strong effect of the plasma on the synchrotron emission at u R * = v p ^/v v jv B . 

Consequently we found that a necessary condition for the collective plasma effects to be 
observable is v v jvB ^> 1. This ratio was shown to depend only on the ratio of £ e to £b, and 
on the width of the distribution function [eq. (7)]. Hence £ e /£,B ^ 1 is a necessary condition 
for the effects to be observable. 

We applied the above results to plasma parameters {pf e ~ 10 3 — 10 4 , v p jvB — 10 2 — 10 3 ) 
typical to GRB afterglows, and found that for ~ 10~ 6 and £ e close to equipartition, the 
plasma had a decisive effect on the synchrotron emission at radio frequencies, during the 
transition of the fireball dynamics to self-similarity. Two expansion scenarios were consid- 
ered : (i) a blast wave into a uniform-density ISM with n e ~ 1 cm -3 , where we expect a 
plasma collective effect on synchrotron emission from the gas heated by the forward shock 
at ~ 3.5 MHz; (ii) expansion into a wind, where the number density decreases as r~ 2 , and at 
the radius at which the dynamics becomes self-similar is ~ 10 4 cm~ 3 . In this case we predict 
a plasma effect on synchrotron emission from the wind accelerated by the forward shock at 




(32) 




0.2 GHz [see eqs. (26)-(29)]. 



Note that plasma collective effects on synchrotron emission below v R * are expected 
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irrespective of the details of the electron distribution function. This is in contrast to neg- 
ative reabsorption and the maser effect, which require a region in the electron distribution 
function which rises faster than E 2 . The energy distribution of electrons produced by col- 
lisionless shock acceleration is not well known, in particular at low energy. However, we 
have shown that for plausible assumptions on the electron distribution function at low en- 
ergy, synchrotron cooling of the electrons may lead to negative reabsorption. This, however, 
depends on the details of the distribution function near its peak. We have demonstrated 
negative reabsorption for a distribution which transforms smoothly (at 7 em i„) from a E 2 
power-law at low energies to a E~ p {p = 2.4) power-law at high energies, for a certain range 
of parameters (Figure 1 Right). 

One should notice, however, that even if the distribution function permits negative 
reabsorption, coherent emission from plasma accelerated by a reverse shock may be com- 
pletely obscured by the high optical depth of the ambient medium heated by the forward 
shock. Thus, for example, when the fireball expands into a wind, certain electron distribu- 
tion functions may lead to coherent emission from the reverse shock at 1.2 GHz [eq. (29)]. 
Nevertheless, the maser due to the reverse shock is completely suppressed by the ~ 10 3 
optical depth of the forward shock at the GHz frequency range. A similar situation occurs 
also when the fireball expands into a uniform- density ISM. Note, however, that if turbulent 
processes disrupt the shock fronts, so that "bulges" of accelerated ejecta lie in front of the 
forward shell along the line of sight, the plasma effects occurring in the reverse shock may 
become observable. 

Our calculations of coherent emission assumed a homogeneous distribution function in 
the emitting shell. However, we have numerically shown that even if inhomogeneity due to 
different cooling times is taken into account, coherent emission would still be observable : 
the regions close to the shock front, where the electrons have little time to cool, contribute 
only a positive optical depth of order of few, hence do not obscure completely the maser 
emission from regions further away from the shock front. 

We have shown that for typical plasma parameters the normal waves propagating in the 
plasma are circularly polarized at frequencies of order of u R * [see the discussion following eq. 
(9) and §3.3]. Since the difference between the refractive indices of the transversal normal 
modes is smaller by a factor of ~ vb/v <C 1 than 1— n\^{v), and it is the latter which enter the 
equations of reabsorption and determine the amplification of radiation in both polarizations, 
the resulting optical depths corresponding to radiation in the two polarizations differ by a 
factor of a few, at most, if the magnetic field is assumed to be constant throughout the 
emitting shell. This factor is expected to reduce, due to saturation of emitted radiation. If, 
however, the direction of the magnetic field is assumed to fluctuate randomly on length-scales 



smaller than the shell's width, the difference between the two refractive indices also changes 
signs randomly, and the difference in the amplifications of the two polarization is averaged 
out. Consequently, we do not expect the observed radiation to have a high degree of circular 
polarization at relevant frequencies. 

A detection of a strong effect in the radio-band at early stages of the afterglow will 
place a new independent constraint on the value of as it will prove that it must be much 
smaller than 1. Since there is a two orders of magnitude difference between our estimates of 
z/r* in the forward shock in the two expansion scenarios, a flash in radio wave-band during 
early afterglow stages seems more probable when the fireball expands into a wind. In this 
respect, plasma collective effects may serve as an important clue to the GRB environment 
and progenitor type. 

EW is the incumbent of the Beracha foundation career development chair. 
A. Dispersion relations of transversal electromagnetic waves in a relativistic 



Since the rate of binary collisions is much smaller than the frequencies of interest for 
us, the plasma may be adequately described by a collisionless Vlasov equation : 



where the one-particle distribution function / a (x Q , p a , t) is defined as the fractional density 
of particles of a single type in phase space, and a denotes a single species of particles, i.e. 
electrons or protons. 



Let f a0 be an equilibrium distribution function, and let f al (x, p, t) be a small pertur- 
bation to it, so that f a = f a o + f a i. It is assumed that the equilibrium distribution function 
carries no net charge or current distributions. Since we assume that \f a i\ f a o, equa- 
tion (Al) can be linearized, and we obtain the (Fourier-transformed) perturbation to the 
equilibrium distribution function 



plasma 




(Al) 



A.l. Transversal EM waves in a field-free plasma 
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where we have used the assumption that f a o is isotropic. Here u 
electromagnetic fields satisfy Maxwell's equations 

ldB 

~c~dt 
IdE An. 



2-ku. The transversal 



V x B 



c dt c ^ ' 



(A3) 



where the net current is given by j = ^2 a n a q a I (p/l m a)faid 3 p, the average number density 
of the particles of the species a is n a , and 7 = ^p 2 + m 2 a c 2 /m a c is the particle's Lorentz 
factor. We choose a system of axes such that the radiation propagates with a wave vector 
k = (0, 0, k) along the z axis, whence the transversal electric field is confined to the xy plane : 
E = (E±i, Ej_2, 0), and the particle momentum vector is p = (p_i_i,p±2,P\\)- Equation (A3) 
leads to a dispersion relation for the transversal waves : 



Arc r/u • Ei = eijEj with i, j = 1,2 



where is the dielectric tensor 



en — C22 



^12 — —£21 
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pa J 



d 3 p- 



dfao P±i m a c 



uj — kv») dp 



P \Jp 2 + m 2 a c 2 







(A4) 

(A5a) 
(A5b) 



Here uj^ = (Annaq^/ma) 1 ^ 2 is the non-relativistic plasma frequency of the species a. The 
refractive index n(u>) = kc/u for the transversal waves is the square root of the doubly- 
degenerate eigenvalue of equation (A4). 

If the electron and proton distribution functions can be approximated by 5— functions 
fao = CL/4np 2 a0 )5(p - p a0 ), equation (A5a) becomes 

,NR\2 



Cll — ^22 



(. ,NH\ 
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+ 2x a ( 1 — 

7a0 



(A6) 



x a + 1 

with x a = u/kv a o, and v a o = Pao/lao^a- It can be shown that in the frequency range of 
interest, and for relevant values of 7 CjP , the second term in the curly brackets of equation 
(A6) dominates the first term for both the electrons and the protons. This may be used to 
approximate the refractive index by : 
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(A7a) 
(A7b) 
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Both observations and theory imply that in afterglow plasmas the electrons are not 
distributed mono-energetically, but rather have a distribution function which extends to high 
energies as a power-law, with a spectral index p > 2. However, we have shown numerically 
that in the frequency range of interest for us, where |1 — (kc/u)\ ^> I/tI, and for relevant 
values of electron and proton Lorentz factors, |1 — upl\ differs from |1 — n$\ by less than one 
part in 10 4 , where n$ and ripi are the refractive indices for the mono-energetic and power-law 
distribution functions, respectively. Consequently, result (A7a) still holds approximately for 
a power-law distribution. Hereafter we shall use 5— function distributions to evaluate the 
deviation of the refractive index from 1. 



A. 2. Transversal EM waves in a weakly magnetized plasma 

Let f a o be an isotropic equilibrium distribution function, which describes a plasma in 
an external uniform magnetic field Bo- We introduce a small perturbation to f a o. Since we 
are interested in frequencies tv ^> ujb, we assume that the external magnetic field is small, 
so that the perturbation is made up of two contributions, one (which we denote by f ai ) 
which is independent of the magnetic field, and the other (denoted by / a2 ) which is linear 
in the magnetic field; thus f a0 3> \f a i\ 3> |/ a 2|- We now linearize the Vlasov equation [see 
eq. (Al)], and neglect (df a2 /dp) with respect to (df al /dp). Substituting equation (A2) for 
fai, we obtain (after some algebra) 



q a (dfgo/dp) 

Ja2 = IP X B Q )- 

^m a c p 



E (E-p)k 



+ 



(A8) 



.(w-k-p/7m a ) 7m Q (u - k • p/jm a ) 

Hence f a 2 is indeed linear in the external magnetic field, consistent with our assumption. 

As already mentioned, we are interested in frequencies uo ^> ujb- For these frequen- 
cies the effect of the magnetized plasma on the propagating radiation is small, and the 
electric field can be assumed to be approximately transverse, i.e. Elk, Our previ- 
ous choice of axes is convenient for the analysis; without any loss of generality, we choose 
B = (0, B sin 0, B cos 0) . Following the steps outlined in equations (A3) - (A5a) , we obtain 
the dielectric tensor for quasi-transversal EM waves in a weakly magnetized plasma : 



ei2 = -e 2 i = % 7 J u Ba COS 



"pa*) 2 , .nr „„„ A f d s p 1 df a0 p 2 ±l mlc 2 

J {uo — kv\\) 2 dp p p 2 + m^c 2 ' 



with ujga = \q a \Bo/m a c the non-relativistic gyration frequency. The expression for en = e 2 2 
are given in equation (A5a). For mono-energetic electron and proton distributions, equation 



(A9) becomes 



E 



,NR\2, ,NR 



Upa ) ^Ba COS</> 2 

f J2 = 1 = '> ■ ^ ^T*« 



a 



aO 



2 1 10 
2 K~ 



where x Q has the same definition as in equation (A6) 
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Fig. 1. — Emitted intensity of circularly polarized synchrotron radiation in plasma. Dashed 
lines represent the "vacuum" values of I v , and solid lines represent these quantities when 
plasma effects are considered. The vertical dash-dotted lines at u R * separate the region where 
plasma effects are negligible [y > v R *) from the region where these effects are important 
[y < v R *). In order to demonstrate the plasma collective effects given different plasma 
conditions, two parameter sets were used for the calculations presented here. 
Left : Typical plasma conditions behind a forward shock propagating into a uniform density 
ISM (y v ~ 2 x IO 5 Hz in the observer's frame, 7 em m = 4 x f0 4 , n' e = IO 3 and ^cooi/lemm — 
2 x f 3 ). Although there is no negative self-absorption, a ~ 1.5 order of magnitude increase in 
the emitted intensity (compared to the emission in vacuum) is apparent at v < u R * ~ 10 3 z/ p 
(see §4.2.3). 

Right : Typical plasma conditions behind a forward shock propagating into a wind {y v ~ 



10 7 Hz as measured in the observer's frame, 7 e 



10 4 , n' e = 10 7 and 7 coo «/7 emi „ ~ 10). 



The divergence of the emission at v < v R * ~ 10z/ p is the signature of the maser effect. 



